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Homogeneous charge compression ignition (HCCI) engine uses a relatively new mode of combustion 
technology. In principle, there is no spark plug or injector to assist the combustion, and the combustion 
auto-ignites in multiple spots once the mixture has reached its chemical activation energy. The 


challenges in developing HCCI engines are the difficulties in: controlling the auto-ignition of the mixture 
and the heat release rate at high load operations, achieving a cold start, meeting emission standards and 
controlling knock. At low engine speeds, early auto-ignition can occur, possibly leading to knocking, 
Keywords: while late auto-ignition at high engine speeds will make HCCI susceptible to misfire. Hydrogen greatly 
HCCI reduces emissions levels but with reduced power. However, when hydrogen is combined with diesel in 
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n dual-fuel mode, low NOx, CO and particulate matter (PM) emissions levels can be achieved, and engine 
Natural gas efficiency can be increased by 13-16%. Numerical methods are commonly used to predict HCCI engines’ 
Diesel performance (i.e. emissions levels, brake thermal efficiency and combustion phasing), which is cost- 
Numerical effective compared to solely relying on experimentation. The multi-zone method promises better 
Chemical kinetic simulation results compared to the single-zone model by combining detailed chemical kinetics with 
simplified 3D modeling so that turbulence and inhomogeneity in the mixture are considered; good 
agreement between simulations and experiments have been achieved. Specific strategies used in the 
experimental method (e.g. fuel additives, inlet air heating, inlet air pressurizing, exhaust gas recirculation 
(EGR) and injection methods), and numerical method (e.g. single-zone and multi-zone models, mixing 
model, turbulence model and multi-dimensional model), and other issues associated with HCCI engines 

are discussed in this paper. 
© 2014 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Internal combustion (IC) engines are widely used in numerous 
applications throughout the world. A new mode of combustion is 
being sought in order to reduce the emissions levels from these 
engines: homogeneous charge compression ignition (HCCI) engine 
technology is a potential candidate. The HCCI technique is the 
process by which a homogeneous mixture of air and fuel is 
compressed until auto-ignition occurs near the end of the com- 
pression stroke, followed by a combustion process that is signifi- 
cantly faster than either Compression Ignition (CI) or Spark 
Ignition (SI) combustion |1]. HCCI technology claimed to improve 
the engine thermal efficiency while maintaining low emissions 
and can be implemented by modifying either SI or CI engines 
using any fuel or combination of fuels [2,3]. The air/fuel mixture 
quality in HCCI engines is normally lean, it auto-ignites in multiple 
locations and is then burned volumetrically without discernible 
flame propagation [4]. Combustion takes place when the homo- 
geneous fuel mixture has reached the chemical activation energy 
and is fully controlled by chemical kinetics [5] rather than spark or 
injection timing. 

Since the mixture is lean and it is fully controlled by chemical 
kinetics, there are new challenges in developing HCCI engines as it 
is difficult to control the auto-ignition of the mixture and the heat 
release rate at high load operation, achieve cold start, meet 
emission standards and control knock [6,7]. The advantages of 
using HCCI technology in IC engines are: (1) high efficiency 
relative to SI engines - approaching the efficiency of CI engines 
due to the ability of these engines to high compression ratio (CR) 
and fast combustion [8,9]; (2) the ability to operate on a wide 
range of fuels [9-11]; and (3) the ability to be used in any engine 
configuration: automobile engines, stationary engines, heavy duty 
engines or small sized engines [2,12,13]. On the other hand, HCCI 
engines have some disadvantages such as high levels of unburned 
hydrocarbons (UHC) and carbon monoxide (CO) [6,14,15] as well 
as knocking under certain operating conditions [6,14,16]. 

Emissions regulations are becoming more stringent and NO, 
and soot emissions levels in HCCI engines have been greatly 
reduced without sacrificing efficiency, which is close to that of CI 
engines [14]. However, knocking is still the major issue because of 
its sudden onset. Knocking is due to premature combustion where 
the ignition takes place before the piston reaches top dead center 
(TDC) and it reduces engine’s reliability due to high vibration 
effects. 

The performance of an HCCI engine is strongly dependent on 
the fuel type, and this affects the emissions levels as well. Since 
the emissions levels become one of the factors driving engine 


technology today, HCCI development has moved to a next level. 
Due to the importance of HCCI technology, which potentially can 
replace the conventional SI and CI engines, there is a need to 
report the recent development of HCCI engines. This paper 
discusses the current issues for this technology. It consists of six 
sections, where Section 2 contains the performance comparison 
between natural gas and hydrogen in HCCI mode. Section 3 
presents the methods to create homogeneous mixtures, with 
ignition control of HCCI engines in Section 4. Section 5 will discuss 
relevant numerical studies, concluding with Section 6. 


2. Performance comparison 
2.1. State of the art current internal combustion engines 


Producing homogenous mixtures to achieve near-complete 
combustion is a common goal that will lead to the development 
of low polluting engines. Some technologies include Fuel-Stratified 
Injection (FSI), Turbo-Stratified Injection (TSI) and HCCI. All these 
technologies (FSI, TSI, and HCCI) are used to improve the combus- 
tion efficiencies by introducing a homogeneous mixture inside the 
combustion chamber. 

FSI and TSI are commercially used by the Volkswagen Aktien- 
gesellschaft (AG), which consists of other child companies: Audi, 
Skoda, Seat, Bugatti, Lamborghini, Bentley and Scania | 17]. FSI uses 
direct injection fuel with high injection pressure, where the 
evaporating fuel has a significant cooling effect on the cylinder 
charge [18]. This effect helps in reducing the knock and therefore 
higher compression ratio can be used. The air enters the combus- 
tion chamber at a certain angle by using a moveable flap inside the 
intake manifold, while fuel is directly injected during the intake 
stroke. The fuel injector is located close to the inlet valve in the 
cylinder head. With the help of the piston crown design, the air 
will experience tumbling effects inside the chamber. This in turn 
will help the fuel mix with air homogeneously. 

TSI engines, on the other hand, implement a high intake 
pressure (using multipoint injectors) on an FSI engine. This allows 
the fuel to mix homogeneously with the air during the compres- 
sion stroke. Furthermore, pressurizing the intake air will assist the 
combustion and therefore produce better efficiency, allowing 
smaller engines to be built with power and torque similar to that 
of bigger engines. The first engine uses direct injection technology 
to stratified charge engine was the Texaco combustion process 
[19], as reported by [20]. 

HCCI engines can be considered as new technology even 
though the research was initially undertaken by Onishi et al. in 
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1979, as reported in [21]. Researchers worldwide are investigating 
HCCI engines as this technology has not yet been sufficiently 
developed and commercially available. They can be used in either 
SI or CI engine configurations with a high CR. HCCI engines operate 
without using diesel injectors or spark plugs and high efficiency 
can be achieved with low NOx and soot emissions. General Motors 
(GM) corporation has unveiled a prototype car with a gasoline 
HCCI engine, which could cut fuel consumption by 15% [22]. The 
engine is able to virtually eliminate NO, emissions and lower 
throttling losses, which improves fuel economy. 


2.2. Fuels used in HCCI engines 


HCCI engines can operate using any fuel as long as the fuel can 
be vaporized and mixed with the air before ignition [10]. Since 
HCCI engines are fully controlled by chemical kinetics, it is 
important to evaluate the fuel’s auto-ignition point to produce 
smooth engine operation: no knocking or misfiring. Different fuels 
will have different auto-ignition points. Fig. 1 shows the intake 
temperature required for different fuels to auto-ignite at difference 
compression ratios when operating in HCCI mode | 10]. It is clearly 
seen that methane requires a high intake temperature and high 
compression ratio to auto-ignite. The auto-ignition point decreases 
when the number of carbon atoms in the hydrocarbon increases, 
as shown in Fig. 1, where methane has the fewest possible carbon 
atoms, while iso-octane has the most of those shown. 

The compositions of natural gas vary for different countries as 
shown in Table 1. It is easily adapted for use as a fuel due its wide 
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Fig. 1. Intake temperature required for fuels to operate under HCCI mode with 
varying compression ratios, reproduced from [10]. 


Table 1 
The difference in natural gas composition between some countries [4,23,35-38]. 


Components %Volume 

Australia Greece Sweden US Malaysia 
Methane (CH4) 90.0 98.0 87.58 91.1 92.74 
Ethane (C2H6) 4.0 0.6 6.54 4.7 4.07 
Propane (C3Hs) 17 0.2 3.12 1.7 0.77 
Butane (C4Hio) 0.4 0.2 1.04 1.4 0.14 
Pentane (C5H12) 0.11 0.1 0.17 - - 
Hexane (CgH 14) 0.08 - 0.02 - - 
Heptane (C7Hi6) 0.01 - - - - 
Carbon dioxide (CO2) 27 0.1 0.31 0.5 1.83 
Nitrogen (N2) 1.0 0.8 122, 0.6 0.45 


Table 2 
Natural gas emissions levels compared to oil and coal (Pounds per billion Btu of 
energy) [24]. 


Pollutant Natural gas Oil Coal 
Carbon dioxide 117,000 164,000 208,000 
Carbon monoxide 40 33 208 
Nitrogen oxides 92 448 457 
Sulfur dioxide 0.6 1122 2591 
Particulates 7.0 84 2744 
Mercury 0.000 0.007 0.016 


availability, is economic and has environmental benefits because it 
produces fewer emissions than crude oil and coal, as illustrated in 
Table 2 [23,24]. It can be seen that crude oil and coal produced 
carbon dioxide at a rate approximately 1.4 to 1.75 times higher 
than natural gas. Natural gas is also more readily available than 
crude oil, with a cost that has been competitive for a long time 
[25]. Natural gas is able to operate as a single fuel in an IC engine 
with low HC and CO emissions, however it suffers less power 
output in HCCI engines [26-28]. Its high auto-ignition point (about 
810 K) gives it a significant advantage over diesel-natural gas 
operation by maintaining the high CR of a diesel engine and 
lowering emissions at the same time [23,29,30]. Duc and Watta- 
navichien [29] claimed that the high octane number of methane 
(about 120) allows engines to operate at a high CR. Results from a 
four-stroke HCCI engine simulation have shown that methane did 
not ignite if the intake temperature was less than 400 K with a CR 
of 15:1, as reported by [31]. This is supported by Fig. 1, where 
methane will only auto-ignite with an intake temperature less 
than 400 K when CR is greater than 18:1. 

Car manufacturers are producing cars powered by fuel-cells, as 
well as engines operated with compressed natural gas (CNG) 
[32,33]. They are purpose-built to reduce emissions and are more 
economical than gasoline and diesel. There are commercialization 
difficulties for fuel-cells because of their complexity to produce 
and safety issues. CNG on the other hand has been successfully 
commercialized with 48% of the global market share, but it has 
difficulties with high capital cost, existing refueling infrastructure 
and competition from other alternative fuels [34]. 

The Indicated Mean Effective Pressure (IMEP) is a measure of the 
engine’s capacity to do work in each stroke, before considering all the 
losses. If IMEP is increased, it can reduce the intake temperature 
needed in an HCCI engine; increasing the CR has the same effect [36]. 
Antunes, Mikalsen and Roskilly [39] have investigated the engine 
performance of an HCCI engine fueled with hydrogen, and they found 
that the IMEP increases when the air intake temperature is decreased. 
Hydrogen requires a lower air intake temperature than natural gas for 
the same CR, as shown in Fig. 1, due to hydrogen having a lower 
density. Therefore, it is important to control the intake temperature 
and the operating conditions of an HCCI engine to achieve an 
optimized combustion point and the methods to achieve this will 
be discussed in Sections 3 and 4. 

Hydrogen can operate as a single fuel in an HCCI engine but it is 
often unstable and is prone to generate knocking [40]. It has the 
highest diffusivity of any gas in air, about 3-8 times faster than 
that of natural gas, which leads to fast mixing [30] and the intake 
charge can be considered homogeneous when premixed with air 
[40]. The net heating value for hydrogen is almost three times 
higher than diesel (119.93 MJ/kg compared to 42.0035 MJ/kg) with 
a high self-ignition temperature to initiate combustion (858 K) 
[41]. Since hydrogen cannot be stably used as a single fuel in an 
HCCI engine, most researchers use it as an additive: either to 
control the ignition timing or to increase the engine performance 
[1,26,30,42,43]. 
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Iso-octane is used as a surrogate fuel for gasoline in HCCI 
engine experiments while n-heptane is used for diesel [44-46]. 
Alcohol-derived fuels, as shown in Fig. 1, are not widely used due 
to their complexity to produce, leading to high manufacturing 
costs. Most of the alcohol-derived fuels like biofuels are still 
subject to worldwide investigations as they create clogging issues 
in the engine. Even though biofuels have their own challenge, they 
have received increasing attention from researchers as a renew- 
able source, as reported by [47]. 


2.3. Natural gas and hydrogen operating conditions 


The combination of Natural gas or hydrogen with diesel 
reported to yield low emissions and to some extent increase the 
engine efficiency, either in HCCI or Cl combustion mode [41,48- 
50]. Diesel alone is not suitable for HCCI engines due to its low 
volatility and high propensity to auto-ignite while natural gas has 
a high resistance to auto-ignition as reported by [51]. 

Combinations of high octane number fuels (such as natural gas and 
hydrogen) with high cetane number fuels (such as diesel) are able to 
increase the engine durability, and under certain operating conditions 
reduce emissions levels such as soot, HC, CO and NO, [40,52-54]. It 
was also reported that these combinations have a high thermal 
efficiency under early injection timing [52,54]. Fuels with a higher 
octane number have better resistance to knocking while fuels with a 
higher cetane number have a shorter ignition delay time, thus 
providing more time for the fuel to complete the combustion. There- 
fore, a combination of both (high cetane number fuels and high octane 
number fuels) provides soft engine run [40], whereby the mixture can 
be operated at high CR and has longer combustion duration. 

Hydrogen has a high octane number (about 130) as well as a 
high lower heating value (LHV) (119.93 MJ/kg) and its combination 
with diesel helps to increase engine efficiency and control the 
auto-ignition point in HCCI engines. Natural gas, on the other 
hand, has a high auto-ignition temperature and it can be used in 
high CR engines. Table 3 compares the physical and chemical 
properties of diesel with natural gas and hydrogen. 

Hydrogen has the highest LHV or Lower Calorific Value (LCV) 
compared to both diesel and natural gas, which means that it releases 
a high amount of energy during combustion and thus produces the 
highest flame speed. The wide range of the flammability limits allows 
a wide range of engine power outputs through changes in the mixture 
equivalence ratio. Flammable mixtures of hydrogen can go from as 
lean as A= 10 to as rich as A=0.14 [50], where J is the air-to-fuel ratio. 


Table 3 
Diesel properties compared to hydrogen and natural gas [30,41,50,55]. 


Properties Diesel Hydrogen Natural gas 
Main component Cy2H23 Hə Methane (CH4) 
Auto-ignition temperature (K) 553 858 923 
Lower heating value (MJ/kg) 42.5 119.93 50 
Density (kg/m?) 833-881 0.08 0.862 
Molecular weight (g/mol) 170 2.016 16.043 
Flammability limits in air 0.7-5 4-75 5-15 

(vol%) (LFL-UFL) 
Flame velocity (m/s) 0.3 2.65-3.25 0.45 
Specific gravity 0.83 0.091 0.55 
Boiling point (K) 453-653 20.2 111.5 
Cetane number 40-60 - - 
Octane number 30 130 120 
CO, emissions (%) 13.4 0 9.5 
Diffusivity in air (cm?/s) = 0.61 0.16 
Min ignition energy (mJ) - 0.02 0.28 


Table 4 
In-cylinder peak pressure and temperature comparison for natural gas and 
hydrogen in various configurations [14,39,56-59]. 


Mode Max pressure (MPa) Max temperature (K) 
HCCI cI HCCI cI 

H2 ~8 ~12 - - 

NG ~7 ~7.5 ~1300 ~1850 
Diesel ~6.1 ~6.6 - ~2300 
NG-+ diesel ~3 ~5.5 ~ 1450 - 
H2+diesel ~7 ~78 - - 

Table 5 


Maximum BTE for hydrogen-diesel fuel compared with hydrogen HCCI, direct 
injection (DI) diesel and dual fuel [57,60]. 


Diesel DI Dual fuel (H2+diesel) Hz+diesel PPCCI H2 HCCI H2 DI 


BTE (%) 27.9 33:9 36.5 45.0 42.8 


2.4. Peak pressure and temperature 


The variations of heat release rate during combustion affect the in- 
cylinder peak pressure and temperature. All these quantities (heat 
release rate, temperature and pressure) depend on the speed of the 
engine, equivalence ratio, load, intake pressure, temperature and 
energy content of the fuel configurations. Higher loads and richer 
mixtures typically produce higher peak pressures. Table 4 shows how 
the in-cylinder peak pressure and temperature vary for HCCI and 
conventional CI modes. These are general data to illustrate which 
mode produces a higher peak pressure. The data for maximum 
temperature is not fully reported but in general, temperature increases 
with pressure. HCCI configurations produce lower peak pressure than 
the conventional CI modes for every fuel configuration, leading to 
significant impacts on emission levels. Hydrogen addition to diesel 
influences the peak pressure generation. The higher the addition, the 
higher the peak pressure, and the addition of hydrogen is able to 
reduce ignition delay as well [40]. It can be concluded that the 
hydrogen fuel configuration yields the highest peak pressure and 
therefore produces better power and efficiency. 


2.5. Brake thermal efficiency 


The brake thermal efficiency (BTE) of an engine is the ratio of brake 
output power to input power and describes the brake power 
produced by an engine with respect to the energy supplied by the 
fuel. A study has been performed [14] to determine the best BTE for 
biogas fuels with an energy ratio ranges from 40% to 57% of biogas and 
an intake temperature of 80 °C, 100 °C and 135 °C. The biogas’ main 
component was methane ( > 60%), which was produced by anaerobic 
fermentation of cellulose biomass materials [29]. It was reported that 
the best energy ratio with diesel in HCCI mode was 51% biogas and the 
optimum efficiency occurred when the intake temperature was 
135 °C. High energy ratios lower the heat release rate and the 
efficiency [14]. However, even when operated at the optimum biogas 
energy ratio, the BTE was no better than diesel running in CI mode. 
Duc and Wattanavichien [29] reported that biogas-diesel running in 
dual-fuel non-HCCI engines has a lower efficiency than diesel single 
fuel in either HCCI or non-HCCI mode. 

Hydrogen, on the other hand, had a higher BTE than pure diesel in 
non-HCCI mode, increasing the BTE by 13-16% [30]. Szwaja and Grab- 
Rogalinski [40] reported that the BTE was increased from 30.3% to 32% 
with an addition of 5% hydrogen. The increase of BTE in hydrogen- 
diesel mode might be due to the uniformity of mixing of hydrogen 
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with air [41]. Hydrogen as a single fuel running in HCCI mode gives a 
BTE of up to 45% [39], showing that hydrogen is able to operate with 
extremely lean mixtures and still maintain a relatively high efficiency 
compared to diesel engines. Table 5 shows that hydrogen in HCCI 
mode yields better results compared to the conventional mode (45% 
vs. 42.8%) and hydrogen with diesel in non-HCCI mode produces a 
higher efficiency than diesel alone [57]. Therefore, it is expected that 
the combination of hydrogen and diesel in HCCI mode will produce a 
higher BTE than diesel in CI mode. This conclusion was reached with 
partially premixed charge compression ignition (PPCCI) configurations 
when reformed exhaust gas recirculation (REGR) was used with 
hydrogen-rich gas (no more than 24%) added to the intake manifold 
[60]. Some researchers use PPCCI instead of HCCI as they differ in their 
injection methods and both of them have the same purpose: to 
achieve a homogeneous mixture. 


2.6. Exhaust gas emissions 


Emissions levels have become a major focus in new engine 
developments nowadays because regulatory bodies, such as those in 
Europe, the United States (US) and Japan, are imposing stringent 
vehicle emissions quality standards [61-63]. Emissions in HCCI 
engines consist of UHC, CO, NOx, and particulates. HCCI engines 
claimed to have low emissions levels on NOx and particulate matters 
[14] and high levels of unburned hydrocarbons (HC) and carbon 
monoxide (CO) [14,64]. However, the emissions levels of an engine 
varies from one engine to another and is dependent on the operating 
conditions of the engine, fuel quality and the engine design [65,66]. 
Thus, the emissions levels from different engines cannot be directly 
compared with another engine because of those factors. 


2.6.1. UHC, CO and CO2 

UHC and CO emissions in HCCI engines are higher than the 
conventional diesel engines as reported in [6,67,68] and are two of 
the challenges that need to be resolved. UHC is the consequence of 
incomplete combustion caused by low combustion temperatures 
[15,69,70] which cause deposition of fuel in boundary layers and 
crevices [58]. The level of UHC is generally specified in total 
hydrocarbon concentration, which is expressed in parts per 
million carbon atoms [71]. The source of UHC is reported from 
the crevice region, cylinder wall with thin layer of oil left when the 
piston moves down and any combustion wall which has a cold 
area [71,72]. Higher concentrations of hydrogen and natural gas in 
diesel engines have the ability to reduce UHC and CO emission 
levels, because the gaseous state of hydrogen and natural gas will 
reduce the wall wetting effect on the cylinder liner [73]. 

The amount of CO, and CO is dependent on the combustion 
efficiency, where the combustion efficiency can be defined as the ratio 
of CO; to the total of fuel carbon present in the exhaust including CO, 
CO, and UHC [74]. CO emissions are controlled primarily by the fuel- 
air equivalence ratio [71,75]. The principal reaction mechanism of CO 
formation is RH—R-— RO, —RCHO—RCO->CO, where R is the 
hydrocarbon radical [71]. The CO oxidation is dominated by the 
reactions of CO with OH radicals as CO+OH = CO2+H, which also 
forms hydrogen radicals. The conversion of CO to CO, occurs when the 
concentration of OH radicals increases during combustion [76]. For the 
HCCI engine with a low combustion temperature, the OH level is 
reduced resulting in an incomplete oxidation of CO [77]. The low 
combustion temperature will reduce the combustion efficiency as a 
result of the lower oxidation activity of the hydrocarbons and the 
lower conversion rate of CO to CO2 [78]. Results of simulations confirm 
that the piston-ring crevice needs to be resolved in order to accurately 
predict UHC and CO emissions [6]. UHC and CO emissions originate in 
the crevices and boundary layer, which are too cold for complete 
consumption [79]. 


2.6.2. NOx and particulate matter (PM) 

The NO, formation is explained by several formation mechan- 
isms: Zel'dovich mechanism, Fenimore mechanism, fuel-bound 
NOx, NO2 mechanism and N20 mechanism. In the Zel’dovich 
mechanism, also known as thermal NO,, the NO, is not formed 
from the fuel because there is no nitrogen component in the fuel. 
NO, is formed in a high temperature reaction, where the nitrogen 
in air dissociates into nitrogen radicals to form NO when reacting 
with oxygen. Some NO is converted to NO2 when further reactions 
occur in the chamber. The thermal NO, is not significant when the 
combustion temperature is below 1800 K [76]. 

Under the Fenimore mechanism, also known as prompt NOx, 
the NOx is promptly formed in laminar premixed flames long 
before the NO, is formed by the thermal mechanism. The Feni- 
more mechanism explains the additional NO, produced over the 
Zel'dovich mechanism in hydrocarbon flames. Prompt NOx is 
important for hydrocarbon fuels in fuel-rich conditions, where 
NO, is formed by rapid reactions of hydrocarbon radicals (CH, CH2, 
Co, CoH and C) with molecular nitrogen. Miller and Bowman [80] 
reported that NO, formed by thermal mechanism is the dominant 
source of NO, only in the equivalence ratio range of @ = 0.8 — 1.0. 
For fuel-air equivalence ratio, @ < 0.8, the temperature is suffi- 
ciently low and NO, was formed by Fenimore mechanism. 

The fuel-bound NO, mechanism is used for coal and coal- 
derived fuels, where nitrogen exists as chemically bound to the 
major fuel. The NO, formation is dependent on the local combus- 
tion temperature, stoichiometric conditions and the level of the 
nitrogen compounds in the fuel-air mixture. The NO2 mechanism 
on the other hand, is based on the chemical kinetic calculations 
near the flame zone, where NO- is formed due to the reaction 
between NO and HOs;. The NO, then reacts with H and O radicals 
to form NOx. The N20 mechanism is also based on the chemical 
kinetic calculations, where the N20 is formed due to the reactions 
of various nitrogen radicals with NO. The N2O will finally react 
with oxygen radicals to form NOx. In short, the NOx formation is 
still under investigation and one cannot claim that all routes have 
been found [76]. 

Generally, most parts of the NO, formation are determined by 
the peak temperature during combustion, where the peak tem- 
perature is dependent on other parameters as well such as 
equivalence ratio, fuel composition and the initial temperature of 
the fuel-air mixture [76]. Tanaka et al. [81] also reported that the 
NO, levels are increased substantially when the equivalence ratio 
is more than 0.33. A higher equivalence ratio causes the maximum 
in-cylinder temperature to exceed 1800 K and produce more NOx. 
In HCCI engines, the NO, and particulate matters (PM) are 
reported to be very low [72,82] by implementing high CR engines. 
NO, can also be reduced using high EGR rates by reducing the local 
temperature and decreasing the oxygen amount in the cylinder 
[70]. However, unburned HC and CO emissions will increase due to 
insufficient oxygen. Bression et al. [70] implemented a high- 
pressure-loop EGR without cooler and Variable Valve Timing 
(VVT) to reduce the NOx, unburned HC and CO by increasing the 
combustion temperature. They reported that internal EGR by using 
VVT is an effective way to reduce the NO, instead of using high- 
pressure EGR. 

Papagiannakis and Hountalas [23] investigated the combustion 
and exhaust emission characteristics of a CI engine fueled with a 
blend of diesel (direct injected) and natural gas (port injected). 
They found that the NO, level for dual fuel engine operation is 
lower than diesel CI engine. The NO, can also be reduced by using 
hydrogen addition discussed in Section 5.3. The concentration of 
NO, between natural gas and hydrogen in diesel HCCI mode is 
different due to different combustion temperatures, because 
hydrogen has a higher temperature and flame speed compared 
to natural gas and diesel [83]. A survey of research papers by 
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Akansu et al. [84] shows that the NOx level is increased when the 
hydrogen content in natural gas-hydrogen mixtures is increased. 
It shows that the combustion temperature and the flame speed of 
hydrogen contribute to a higher level of NO, emissions. 

In biogas—diesel HCCI engines, the NO, level was low when the 
biogas energy was increased [14]. This might be due to a higher 
homogeneity level achieved between air and fuels. Van Blarigan 
[85] in his study reported that the mixture must be homogeneous 
and lean in order to eliminate the production of NOx. Olsson et al. 
{36] stated that the NO, level is low in natural gas HCCI engines 
and when combined with exhaust gas recirculation (EGR), it drops 
further [15]. Even in natural gas—diesel non-HCCI mode, the NO, 
level is lower than diesel conventional CI engines [86]. Hydrogen 
on the other hand produces zero UHC, CO and CO,, due to the 
absence of carbon in the fuel, but still produces NO, [50]. 
Hydrogen operated as a single fuel in CI mode yielded lower 
NO, levels than diesel [57]. Saravanan and Nagarajan [30] in their 
study of hydrogen addition to diesel in non-HCCI mode showed 
that lower NO, emissions were obtained for all load ranges 
compared to diesel in conventional mode. They reported that the 
formation of NO, depends on temperature more than the avail- 
ability of oxygen. Therefore, hydrogen-diesel in HCCI mode results 
in extremely low NO, emissions levels with no significant amount 
of PM [87]. 

Kayes and Hochgreb [88] developed a mechanism model for 
PM and they reported that PM is produced due to liquid fuel 
nucleation from port fuel injection. The PM can also be produced 
from gas-phase nucleation in fuel-rich conditions. The air-to-fuel 
ratio strongly affects PM and burning liquid fuel will increase the 
amount of PM in the exhaust gas emissions [89], hence the use of a 
high pressure injection is advisable so that the fuel is in an 
atomized condition. PM can be reduced significantly when hydro- 
gen is added to a diesel fueled engine [49] and also by heating the 
inlet [88]. When the engine runs in a dual fuel mode between 
diesel and natural gas, no significant reduction in PM was reported 
[90]. Hydrogen addition, on the other hand, shows a good 
advantage in reducing the PM emission levels. 


2.7. Brake mean effective pressure 


The brake mean effective pressure (BMEP) is an effective 
comparison tool to measure engine performance and indicates 
an engine’s capacity to produce power output over the full engine 
speed range. It is also used to compare one engine’s performance 
with another. A high BMEP shows the ability of the engine to 
perform high load operations. One of the HCCI engine’s challenges 
is its limited load range because high load operations tend to 
produce knock [51]. Table 6 shows the BMEP ranges of standard 
engines. In comparison, the BMEP range for passenger cars (Saab 


Table 6 
BMEP range for various engines types [71]. 


Engine type Compression ratio BMEP range (MPa) 
SI engines 

Small (motorcycles) 6-11 0.4-1 
Passenger cars 8-10 0.7-1 
Trucks 7-9 0.65-0.7 
Large gas engines 8-12 0.68-1.2 
Diesel engines 

Passenger cars 17-23 0.5-0.75 
Trucks 16-22 0.6-0.9 
Large trucks 14-20 1.2-1.8 
Locomotive 12-18 0.7-2.3 
Marine engines 10-12 0.9-1.7 


5 cylinders engine) operating in HCCI mode is only between 0 and 
0.36 MPa [91], which is lower than the conventional SI engines. 

Natural gas-diesel in HCCI mode operates in the BMEP range of 
0.25-0.4 MPa [14], which is low compared with the engines in 
Table 6. To increase the BMEP, the concentration of natural gas has 
been increased but if its concentration compared to diesel is too 
high, this will lead to knocking [14]. Fig. 2 shows that the best 
energy ratio for natural gas—diesel HCCI mode is about 51% and 
within the shaded region. For this optimized energy ratio, the 
maximum BMEP is only 0.4 MPa, which is a very limited load 
range and is not suitable for high load engine operations. 

Hydrogen with diesel in HCCI mode is stable up to 0.6 MPa 
BMEP [60,87]. If a supercharger is used in hydrogen-diesel dual- 
fuel mode on a non-HCClI engine, a maximum BMEP of 0.91 MPa 
was reported [92]. Therefore, hydrogen and diesel in HCCI mode 
might be able to be used for high load engine operations. 


2.8. Knocking 


Knocking in SI engines is a phenomenon where the unburned 
mixture in the compressed gas ignites before it is reached by the 
propagating flame front [93]. Knock is physically detected when 
the engine vibrates excessively and a pinging sound can be heard 
outside as a result of the combustion activity. It causes loss of 
power and if not controlled, knocking could lead to severe engine 
damage and shorten its life. Knocking can occur in any reciprocat- 
ing engine. HCCI engines are prone to knock since they are 
controlled by chemical kinetics and there is no fixed mechanism 
to control knock in them. Knocking phenomena limits the load 
range of an HCCI engine: high load operations can easily initiate 
knock, so upper load limits have to be applied [15]. In all engines, 
knocking occurs when the combustion starts before the piston 
reaches TDC, while misfire is when combustion commences after 
TDC. Knocking and misfire are two different behaviors which must 
be avoided in engine operation as both of them can contribute to 
deterioration of engine performance. 

In natural gas-diesel HCCI engines, a misfire occurs when a 
high natural gas flow rate is combined with a low diesel flow rate 
while knocking starts with the opposite configuration [14]. If HCCI 
engines operate on hydrogen-diesel fuels, knocking is expected to 
occur when high amount of hydrogen is added. Knocking will take 
place if the hydrogen content is more than 16% of the energy ratio 
[40] which is confirmed by Guo et al. [94], who reported that it 
was necessary to use hydrogen with mass fraction less than 15% to 
achieve stable combustion. If an HCCI engine runs on natural gas 
with dimethyl ether (DME), the knock limit is at an in-cylinder 
pressure of 9 MPa and a very limited load range is obtained while 
getting unstable operation for high natural gas concentrations. 
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natural gas and diesel. The shaded region shows the best energy ratio for natural 
gas-diesel HCCI engines, reproduced from [14]. 
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This was discovered by using Computational Fluid Dynamics (CFD) 
simulations with detailed chemical kinetics [51]. An investigation 
of auto-ignition and combustion of natural gas HCCI engines 
showed that knocking starts to happen when the equivalence 
ratio is less than 0.45 with an intake temperature of 380 K [16]. 
Fig. 3 shows the area of knocking and misfiring of natural gas HCCI 
engines with constant intake pressure. Thus, natural gas HCCI 
engines will not be able to operate at high load conditions due to 
this limitation. 

Knocking occurs when a rapid release of energy in the 
remaining unburned mixture causes a rapid increase in local 
pressure. In hydrogen-diesel engines, the engine will not work 
in a stable condition if 100% hydrogen is used. A combination of 
hydrogen and diethyl ether (DEE) produces severe knock on 100% 
load [41]. As such, higher hydrogen content leads to a higher 
probability of knock onset. Knocking phenomena can be detected 
in the in-cylinder pressure variations, observing a rapid instanta- 
neous local pressure rise. The graph formed depends on the 
knocking frequency: the higher the frequency, the more severe 
the knock [40]. Fig. 4 shows knocking behavior in a hydrogen HCCI 
engine for a single combustion event with a compression ratio of 
17 by monitoring its in-cylinder pressure. It can be seen that the 
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Fig. 3. Knocking and misfiring area for natural gas HCCI engines under constant 
intake pressure, reproduced from [16]. 
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Fig. 4. Knocking phenomena in a hydrogen HCCI engine for a single combustion 
event, reproduced from [40]. 


in-cylinder pressure suddenly increases to the peak before follow- 
ing the normal pressure trend: the peak pressure is higher than 
the normal in-cylinder peak pressure without knocking. The 
maximum increase in pressure for knocking is unpredictable and 
this abnormality usually occurs when there is an incorrect operat- 
ing condition. Flowers et al. [26] have described the HCCI ignition 
by H20- decomposition. It accumulates due to low temperature 
reactions in the compression stroke. Once the in-cylinder tem- 
perature reaches 1050-1100 K, H202 decomposes rapidly into two 
OH radicals forming an enormous amount of OH radicals, which 
will rapidly consume the fuels. Reduction in OH concentration in 
low temperature regions leads to incomplete combustion, which 
delays the high-temperature oxidation [95]. However, there is no 
concrete mechanism in the literature describing the knocking 
phenomenon in HCCI engines and the methods to control it. In 
SI engines, it is quite easy to control knock. For example, car 
manufacturers usually install knock and oxygen sensors in the 
engine. If the mixture is prone to generate knocking, the sensors 
will automatically adjust the spark timing accordingly. However, in 
HCCI engines, there is no spark to control such situations and there 
is a strong dependence on the right mixtures and conditions. 


3. Diesel injection methods for homogeneous mixtures 


Vaporized diesel fuels are needed to initiate combustion in the 
conventional diesel CI engines, therefore a high diesel injection 
pressure is required. A higher injection pressure leads to better 
fuel vaporization and forms many small droplets of fuel leading to 
efficient fuel-air diffusion processes. Compression of air generates 
high in-cylinder temperatures, which is important to achieve the 
auto-ignition point of the diesel fuel. CI engines using common- 
rail injection systems employ injection pressures of at least 
35 MPa to ensure the injected fuel is vaporized before the 
combustion takes place. However, in HCCI combustion using 
diesel, it is quite challenging to form a homogeneous fuel-air 
mixture before the initiation of combustion due to diesel’s low 
volatility. Therefore, there are three methods to achieve homo- 
genous mixtures in HCCI engines [21]: 


1. Port injection. 
2. Early injection. 
3. Late injection. 


3.1. Port injection 


Port injection is the method of locating the diesel injector in 
the intake manifold, where the injection pressure is not as high as 
a direct injection system. Ganesh and Nagarajan [67] have studied 
the HCCI combustion of diesel with port injection, which they 
called “external mixture formation”. It is used to form a homo- 
geneous diesel-air mixture with injection pressures of no more 
than 0.6 MPa. 

As the injection pressure is not high, the diesel fuel will not 
vaporize when entering the inlet manifold, creating undesirable 
inhomogeneity problem. To overcome this, the fuel vaporizer is to 
be mounted in the air intake manifold system so that the diesel is 
supplied in vapour form [67,69]. However, there are drawbacks 
when using port injection systems to form a homogeneous 
mixture: the soot is quite significant and NO, is not as low as SI 
engines [21]. This might be due to the poor vaporization of the 
diesel, creating inhomogeneous mixtures. Another drawback is 
that the compression ratio has to be reduced to overcome knock, 
which leads to a reduction in engine efficiency. At the same time, it 
is difficult to control the start of combustion. 
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3.2. Early injection 


Early injection uses the direct injection method, but with 
injection timing earlier in the compression stroke; there is no 
modification of the injection pressure from CI engines. Garcia et al. 
[58] have studied the performance of diesel engines operating in 
HCCI mode with early injection by varying the injection timing 
between 17° and 144° crank angle (CA) before top dead center 
(BTDC). It was found that CO and HC were higher in HCCI mode 
compared to the conventional CI mode due to fuel impingement 
on the cylinder walls [58]. Another study of the spray structure of 
common-rail type high-pressure injectors in HCCI engines 
obtained similar results [96]. The spray structure was simulated 
for early injection timing between 240° and 340° CA BTDC with an 
increment of 20° CA with injection pressures of 80 MPa and 
100 MPa. It was found that if the fuel was injected earlier than 
280° BTDC, the fuel would impinge on the cylinder wall and create 
a rich mixture, leading to incomplete combustion due to wall and 
crevice effects. However, if the fuel was injected at 300° CA BTDC, 
the fuel would impinge on the outer surface of the piston area. Any 
injection earlier than 320° CA BTDC, results in impingement on the 
piston bowl. This process is shown in Fig. 5. 

It is obvious that early injection systems lead to fuel impinge- 
ment problems on the cylinder wall and crevice areas. There are 
few methods to address this issue: using multiple injectors, multi- 
pulse fuel injection or narrow spray angle injectors. Multiple 
injectors will cause additional costs while multi-pulse fuel injec- 
tion requires very precise control over injection timing by mana- 
ging the injection pulse and its period [21]. Therefore it is hard to 
control the injection for different engine loads. Narrow spray 
angles might be an easier and cheaper way, which can be achieved 
by replacing the conventional injector (spray angle of about 156°) 
with a narrow spray angle injector (of about 60°). Kim and Lee [97] 
have studied the effect of narrow spray angle injectors in an HCCI 
diesel engine with an injection pressure of 100 MPa. They found 
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Fig. 5. Spray structure simulation result for various initial conditions, reproduced 
from [96]. The dark curved areas at the center represent the piston crown surface, 
the black circle line is the piston crown edges, the outer circle is the cylinder wall 
which in turn resembles piston edge and the five radial lines represent fuel being 
injected from the center. Spray structure behavior does not change much with 
initial conditions while it does with injection timing. 


Fig. 6. Schematic diagram of the narrow spray angle injector, reproduced from [97]. 


that the narrow spray angle injector in HCCI engines using early 
injection timing is effective in maintaining a high IMEP [97]. The 
schematic diagram for narrow spray angle injectors is shown in 
Fig. 6. By using narrow spray angle injectors, the fuel will hit the 
piston bowl instead of the cylinder wall. However, the obvious 
drawback of using early injection is the difficulty in controlling the 
start of combustion. Therefore, one has to pay extra attention to 
the ignition delay effects. 


3.3. Late injection 


Late injection still uses a direct injection approach in order to 
minimize modification to the injection system. However, a high- 
injection pressure system is used with retarded injection timing at the 
end of the compression stroke or at the early power stroke. Conven- 
tional diesel ignition timing varies from 23° BTDC to 5° after top dead 
center (ATDC) depending on the equivalence ratio [98]. For late 
injection systems, typically the ignition timing is retarded up to 7° 
ATDC [21]. Late injection systems require a long ignition delay and 
rapid mixing rate [21,67]. The former can be achieved by implement- 
ing EGR (as discussed in Section 4.4), while the latter is dependent on 
the geometry of the combustion chamber. The Nissan Motor Company 
has developed a successful late injection system for HCCI engines by 
using Modulated Kinetics (MK) combustion systems [21]. They 
implemented high swirl effects in a toroidal combustion-bowl to 
achieve rapid mixing with a bigger piston bowl diameter. A higher 
piston bowl diameter is better as it can minimize the potential of fuel 
impingement on the piston bowl and the cylinder walls; Fig. 7 shows 
an overview of a piston bowl with toroidal geometry. This late 
injection system is able to control the start of combustion better than 
the first two methods. 
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Fig. 7. Piston bowl with toroidal geometry, reproduced from [99]. 


4. Ignition control in HCCI engines 


Ignition control is one of the challenges in developing HCCI 
engines. The challenges include [6,62]: (1) auto-ignition control; 
(2) limiting the heat release rate at high load operations; (3) meet- 
ing emissions standards; (4) providing smooth engine operation 
by achieving cold starts (‘startability’ of the engine); and (5) limited 
load range. Various studies have suggested that ignition can be 
controlled by using promoters or additives, blending of low cetane 
number fuels with high cetane number fuels [100], pre-heating of 
the intake air [15], pressurizing the intake air [101], hydrogen 
addition and varying the amount of exhaust gas recirculation 
(EGR) [26] through early closure of the exhaust valve. 

Managing the ignition delay is one of the effective ways to 
control the start of combustion. Ignition delay is the time lag 
between the start of injection and the start of combustion. This is 
difficult to control because the combustion in HCCI engines is fully 
controlled by chemical kinetics. Too short an ignition delay leads 
to knocking while too long an ignition delay leads to misfiring. The 
ignition delay is strongly dependent on the gas temperature: an 
increase in temperature accelerates the chemical reactions, 
thereby helping to reduce the ignition delay [57,102]. Tanaka 
et al. [103] studied two-stage ignition in HCCI combustion, 
suggesting that the ignition delay can be controlled by the fuel- 
air ratio, initial temperature and additive dosages. 


4.1. Pre-heat intake air 


As discussed in the previous section (in Section 2.2), using 
natural gas and hydrogen in HCCI require a high intake tempera- 
ture. By pre-heating the intake air, ignition delay is reduced and 
thus the ignition timing can be controlled. As shown in Fig. 1, 
methane and natural gas require high temperatures to auto-ignite 
and methane did not ignite for an intake temperature of 400 K at 
low CR [100]. Therefore, it is imperative to pre-heat the intake air 
to make the fuels combust smoothly. 

Antunes, Mikalsen and Roskilly [39] stated that heating the inlet 
air is the most useful method to control the ignition timing. Questions 
may arise regarding the practicality of including electric heaters [16] in 
an engine bay just for this purpose: the heater causes the operation 
and maintenance costs to increase and contributes to extra engine 
weight. Installing a heater is an option used by most researchers as it 
is the easiest way to get intake air heated to some specific 


temperature. Exhaust gas recirculation (EGR) could be another option 
to reduce the need for a high intake temperature [15,43,69]. 


4.2. Pressurized intake air 


Turbochargers and superchargers are commonly used in real 
engine applications because they can be applied to any internal 
combustion engine. The operational concepts of these two devices 
are the same: to provide a high intake pressure into the combus- 
tion chamber, increase the charge density and thereby increase the 
engine performance. Some studies [101,104,105] show that the 
start of combustion (SOC) is advanced if the intake pressure is 
increased by 0.1 MPa. This indicates that pressurized intake air is 
able to improve the auto-ignition of the fuel. However, these 
situations also depend on the type of fuel used, and in this case 
they used primary reference fuels and gasoline. By increasing the 
intake pressure, it was possible to get the auto-ignition to start at 
15° before top dead center (BTDC). On the other hand, super- 
charging (pressurizing intake air) is able to increase engine 
efficiency [7]. A supercharged hydrogen-diesel engine, but in 
non-HCCI mode, was able to maintain high thermal efficiencies 
and it was possible to use more than 90% hydrogen energy 
substitution for the diesel [60]. Another study by Guo et al. [94] 
shows that hydrogen-diesel in HCCI mode with pressurized intake 
air (150 kPa) is able to improve the atomization process, and 
therefore improve the combustion efficiency. CO, HC and NO, 
emissions were also decreased in this case. 


4.3. Hydrogen addition 


Hydrogen is one of the promising renewable fuels because it is 
naturally available on the earth and can be produced from various 
resources such as fossil energy and biomass [106,107]. Hydrogen 
can be used as a sole fuel in an SI engine but cannot be used in a CI 
engine due to its higher ignition point compared to diesel 
[108,109]. The used of hydrogen in internal combustion engines 
requires extra care due to safety reasons. According to Najjar [106], 
there are three hazards associated with the use of hydrogen, 
which are physiological (frostbite and suffocation), physical effect 
(embrittlement and component failures) and chemical (burning or 
explosion). Proper hydrogen installation such as high pressure 
vessel is required to avoid physiological effect. The ability of 
hydrogen to permeate through steel may cause an embrittlement 
effect and lead to component failures. Fayaz et al. [110] reported 
that to avoid the embrittlement effect on some metals, the 
following alloys can be used in a hydrogen engine application: 
brass and copper alloys, aluminum and its alloys and copper 
beryllium. The crank case design of the engine should also be 
taken into consideration so that it has good ventilation. The blow- 
by effect causes unburned hydrogen to enter the crankcase and at 
certain concentrations may lead to combustion. The hydrogen 
injection method also plays an important role in avoiding unde- 
sired explosions. Direct injection is a preferred method over port 
injection to avoid the backfire effect [110]. Backfire is combustion 
occurring during the intake stroke due to hot spots and could also 
occur in the intake manifold [50]. It is of the utmost importance to 
avoid any undesired combustion behavior for safety reasons. 

Hydrogen can be produced in many ways including natural gas 
conversion, coal gasification, electrolysis, biomass gasification and 
photolytic processes |110]. However, the conversion cost of hydro- 
gen using current technology is still very expensive [106] and it is 
not practical to have a separate hydrogen vessel for a vehicle with 
dual fuel technology because the combined fuel cost becomes 
higher. Thus, an on-board hydrogen generator is a possible solu- 
tion to overcome this high hydrogen conversion cost and also the 
problems associated with storage and handling [111]. Tsolakis and 
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Megaritis [87] used a fuel reforming reactor in the diesel exhaust 
system to extract hydrogen. Bromberg [112] used a plasmatron 
fuel converter to convert a gas containing hydrogen, while Bari 
and Esmaeil [113] used water electrolysis to generate hydrogen. 

Hydrogen addition in a CI engine is able to reduce UHC, CO and 
smoke emissions, however, NO, emission is increased due to the 
high combustion temperature [49,107,109-111,114,115]. To reduce 
the NO,, intake charge dilution with nitrogen can be used [111], or 
EGR [110]. With nitrogen dilution, the NO, level is reduced 
substantially, however, smoke, CO and fuel consumption are 
increased. EGR, on the other hand, causes increase in PM, UHC, 
CO and also engine wear. Thus, Miyamoto et al. [116] recom- 
mended using large amount of EGR and late diesel-fuel injection 
timing for diesel engine with hydrogen addition. They reported 
that the NO, emission is low without the increase in UHC due to 
low combustion temperature. 

Hydrogen addition to the HCCI engine is one of the effective 
ways to reduce the ignition delay due to its high diffusivity in air, 
which causes rapid mixing between fuels and oxidizers. Hydrogen 
addition in a natural gas mixture is able to increase the in-cylinder 
peak pressure, reduce ignition delay time and ignition tempera- 
ture, and increase indicated power [1]. It also allows the extension 
of the lean limit of the natural gas mixture, without entering the 
lean misfire region, while achieving extremely low emissions [50]. 

Hydrogen addition in ultra-low sulfur diesel (ULSD) promotes 
partially-premixed compression ignition and results in improved 
performance and reduction in emissions [60]. In HCCI engines, 
hydrogen is able to form a homogeneous mixture due to its rapid 
mixing behavior, and thus increase the efficiency of the engine 
{41]. Furthermore, hydrogen can be produced from the exhaust 
gases of the engine itself using a reformer, which is called an “on- 
board hydrogen producer” [15,87]. As the amount of hydrogen is 
increased, the auto-ignition delay time reduces accordingly while 
as the in-cylinder peak pressure is increased, the ignition tem- 
perature reduces and indicated power increases [1]. The addition 
of hydrogen is not costly because it uses a lower-pressure fuel- 
injection system [79]. 

Furthermore, the addition of hydrogen increases the engine 
efficiency by a significant margin, about 13-20% [30,41]. By using a 
catalytic reforming aid in HCCI, the addition of hydrogen in natural 
gas HCCI engines helps in decreasing the need for high intake 
temperatures and also is a means of extending the lower limit of 
HCCI operations [43]. It can be seen from Fig. 8 that the in-cylinder 
peak temperature increases and the ignition delay time reduces as 
the amount of hydrogen increases. 

Szwaja and Grab-Logarinski [40] studied hydrogen addition (in 
HCCI mode) with diesel in a CI engine and found that the addition 
of hydrogen in small amounts (e.g. about 5% in energy ratio) was 
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Fig. 8. Hydrogen effect on the in-cylinder peak temperature and the ignition delay [1]. 


able to reduce the ignition delay and improve engine performance. 
They also concluded that the addition of hydrogen to diesel should 
not be more than 15% in energy ratio to avoid severe knock. 


4.4. Exhaust gas recirculation (EGR) 


EGR is a system where some portions of the exhaust gas are 
returned to the combustion chamber, controlled by using a valve. 
Exhaust gas contains less oxygen (O2) and more CO: combustion 
with reduced O; will lower the combustion temperature, and thus 
reducing the NO, emissions. It is also used to reduce the large heat 
difference between peak pressure and intake pressure. This means 
that the intake temperature of using EGR is high and certainly 
helps in auto-ignition of the mixture. Saravanan and Nagarajan 
[30] reported that EGR is able to increase the BTE of the engine, 
while reducing the NO, formation. Furthermore, EGR improved 
auto-ignition of the engine and reduced the in-cylinder peak 
pressure. EGR used in a diesel HCCI engine showed an improve- 
ment of 1.1% engine efficiency, advanced the auto-ignition by 10° 
and reduced the heat release rate by 11 J/deg CA compared to 
diesel fuel [69]. Ganesh and Nagarajan [67] found that by using 
EGR on diesel HCCI engines, they were able to reduce the heat 
release rate with a lower combustion temperature. They could 
control ignition delay as well as reduce the in-cylinder peak 
pressure. However, using EGR leads to a power loss of 20% when 
used with diesel fuels [58]. 


4.5. Addition of dimethyl ether (DME), hydrogen peroxide (H202), 
formaldehyde or n-butane 


Additives do not seem to lead to any major practical improve- 
ments, due to the addition of extra chemical compounds, which 
leads to an increase in operational cost. However, with the help of 
additives, the start of combustion of DME in a natural gas/DME 
mixture is improved by reducing the ignition delay time. It was 
found that methane-fueled HCCI engines with an intake tempera- 
ture of 400K ignited only if small amounts of additives were 
added [100]. This shows the advantage of using additives, which 
can control the auto-ignition timing, but with extra drawbacks. 


5. Numerical study of HCCI engines 


Simulations are undertaken to reduce research costs while 
maintaining good productivity because of their cost efficiency 
compared to experimentation alone. By using numerical methods, 
one can optimize engine parameters before conducting experi- 
ments and obtain optimized parameters within a short time and 
low cost. Many researchers use KIVA-3V CFD software in combina- 
tion with detailed chemical kinetics solutions using CHEMKIN 
from Sandia National Laboratories [104,117-120]. Good agreement 
between simulations and experiments has been achieved. 


5.1. Chemical kinetics 


In a CI engine, the fuel is direct-injected in the chamber with high 
injection pressure when the piston is nearly at TDC. Then, the fuel 
ignites rapidly in the hot air environment. In SI engines, the spark plug 
is triggered when the piston is approximately at TDC to initiate the 
combustion. HCCI engines, on the other hand, have no mechanism to 
control the ignition timing and rely solely on chemical kinetics for 
combustion, as discussed in Section 3. The combustion in an HCCI 
engine is triggered when the heat in the chamber has reached the fuel 
activation energy. To numerically investigate the combustion behavior, 
chemical kinetic mechanisms which represent the actual fuels have 
been developed [121-126]. A detailed mechanism such as those 
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Fig. 9. Ignition delay time for n-heptane validated against experiment for different 
pressures [131]. 
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Fig. 10. Summary of ignition burn rate and delay time for different types of fuels [103]. 


developed for diesel fuel consists of a blend of components such as 
n-alkanes, iso-alkanes, cyclo-alkanes, aromatics and others [127]. 
Those components were then combined to represent the real fuels 
[128]. Most studies on developing the chemical kinetic mechanism 
[46,121,122,125,127,129] for HCCI engines use the Arrhenius-like plot 
(such as shown in Fig. 9) to relate the fuel activation energy with 
ignition delay time, where the activation energy is proportional to the 
slope of the graph [130]. 

Small variations in fuel properties affect the HCCI combustion 
by having a different ignition delay time. Fig. 10 shows the ignition 
delay time and burn rate of different fuels [103]. The n-heptane 
has the lowest ignition delay time, while iso-octane has the 
highest with a slower burn rate. This shows that the ignition delay 
time (for multiple ignition stages as in Fig. 10) and burn rates are 
dependent on the fuel structure. Tanaka et al. [103] reported that 
the failure on Toluene auto-ignition is due to the high activation 
energy required to abstract the hydrogen (H) radical from the 
carbon chain. 


Westbrook et al. [127] studied the detailed chemical kinetic 
models for large n-alkanes, which can be used to represent the 
diesel fuel. The aim was to investigate the model that can use less 
computing resources by having an almost identical ignition 
behavior compared to the experiment. They found that the 
ignition delay time is not affected significantly for different carbon 
lengths. Thus, they concluded that a small n-alkane such as 
n-octane can be used as a surrogate fuel for diesel which has the 
same ignition behavior of a much larger n-alkane like n-hexade- 
cane. Many studies have been conducted to develop a chemical 
kinetic mechanism that can represent the actual fuels. This is 
summarized in Table 7, where it shows different chemical kinetic 
mechanisms to represent different fuels. 

It is common to use n-heptane as a surrogate fuel for diesel 
[132-135] because the chemical properties between those two are 
very similar in term of cetane number. The intake temperature 
when n-heptane is used as a fuel in a HCCI engine is not as high as 
methane or natural gas. A study by Guo et al. [132] showed that 
the n-heptane can be easily ignited when the inlet temperature is 
313 K on a low CR engine (CR=10). Methane, on the other hand, is 
reported to have no ignition when the inlet temperature is less 
than 400 K on a high CR engine (CR=15) [31]. Thus, it is important 
to have the right fuel and its chemical kinetic mechanism in order 
to study the combustion behavior of an HCCI engine. 

Curran et al. [121] developed a detailed mechanism for 
n-heptane and validated the result over a wide range of operating 
conditions. They found that the ignition delay is in very good 
agreement with experiments using flow reactors, shock tubes and 
rapid compression machines. Because a detailed mechanism uses 
more computational resources, Seiser et al. [136] and Golovitchev 
[123] developed a reduced mechanism for n-heptane. Patel [137] 
reduced the mechanism by more to obtain 26 species and 52 
reactions. He reported that the CPU time is reduced by 50-70% and 
the validation was completed under both constant-volume reactor 
and HCCI engine conditions. The ignition delay result is similar to 
those of detailed mechanisms. The use of n-heptane as the only 
main component in the mechanism is not sufficient to predict the 
soot emission from CI and HCCI engines. Thus, Wang et al. [138] 
improved the mechanism by blending the n-heptane and toluene 
mechanisms. They reported that the mechanism gives reliable 
soot predictions and combustion phasing under various engine 
conditions. 

An iso-octane mechanism was used as a surrogate fuel for 
gasoline and it can also be blended with n-heptane to represent a 
real fuel [81,122]. A real gasoline fuel consists of thousands of 
hydrocarbon compounds [139,140] and the fuel should have been 
modeled using a combination of a few components. Mehl et al. 
[125] developed a mechanism that represents a commercial grade 
gasoline, which consists of n-heptane, iso-octane, toluene and 
olefins. They reported that the result is in good agreement over a 
wide range of pressures and temperatures relevant to internal 
combustion engine applications. A reduced mechanism was also 
developed by Mehl and Chen et al. [141] and further reduced by 
Lee et al. [124]. A reduced mechanism shows an advantage when 
used in a complex CFD model, which requires more computational 
resources [124]. 

Natural gas mainly consists of methane as the main compo- 
nents, as shown in Tables 1 and 3. A mechanism for methane 
called Gri-Mech [142] has received wide attention, such as 
[31,100,143]. A further development of the methane mechanism 
has ended with the final version being Gri-Mech 3.0. Then, 
reduced mechanisms were developed based on the Gri-Mech 
mechanisms such as those from [144]. A recent development of 
the methane mechanism has been developed by Slavinskaya et al. 
[145]. The mechanism was validated for high- and low-pressure 
conditions for rocket applications. They take into consideration 
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Table 7 
Chemical kinetic mechanisms available for different fuels. 


Fuel Main components Year 
Diesel n-Heptane 1998 [121 
n-Heptane 2000 [136 
n-Heptane 2003 [123 
n-Heptane 2004 [137 
n-Heptane/toluene 2013 [138] 
Gasoline Iso-octane 2001 [139 
Iso-octane 2002 [122 
Iso-octane 2002 [140 
Iso-octane 2003 [81] 
n-Heptane/iso-octane/toluence 2003 [123 
n-Heptane/iso-octane/toluene/olefins 2011 [125 
n-Heptane/iso-octane/toluene/olefins 2011 [141 
n-Heptane/iso-octane/toluene 2011 [124 
Natural gas Methane 1999 [142] 
Methane 2001 [151] 
Methane 2006 [152 
Methane 2007 [153 
Methane 2012 [144 
Methane 2013 [145] 
Hydrogen Hydrogen 1995 [146] 
Hydrogen 2004 [154 
Hydrogen 2012 [155] 
Hydrogen 2012 [147 
Biofuels Ethanol 1999 [156] 
Methanol 2002 [157 
Ethanol 2006 [158 
Methyl butanoate 2007 [159 
Methyl butanoate/n-heptane 2008 [160 
Methyl decanoate/methyl-5- 2010 [148 
decenoate/ 
methyl-9-decenoate/n-heptane 
Methyl decanoate/methyl-9- 2012 [161] 
decenoate/ 
n-heptane 
Butanol 2012 [162] 
Methyl decanoate/methyl-9- 2013 [163] 


decenoate/n-heptane 


the depletion of ozone layer. The result shows that a methane-air 
propelled liquid rocket engine influences the formation of nitrogen 
compounds and thereby depleted the ozone. 

Hydrogen mechanisms, on the other hand, show a matured 
development, where the species and reactions are not substan- 
tially varied. A well-developed hydrogen mechanism in 1995 by 
Marinov et al. [146] has been validated over a wide range of 
operating conditions, including high-pressure combustion. A 
recent development by Burke et al. [147] focused on the discre- 
pancies (in reaction rate parameters) between experiment and 
simulation. They reported the characterization of the non-linear 
mixture for H+0O2(M)=HO2(+M) might be necessary to predict 
high pressure flame speed within an accuracy of 5%. 

Biofuel is considered as an alternative fuel and is produced 
from long-chain fatty acids derived from either vegetable or 
animal oil. The trans-esterification method is used to transform 
these large molecules into esters, where methanol is commonly 
used during the trans-esterification process [148]. The focus of 
using biofuels has received increasing attention [47] because of the 
advantage as alternative fuels. Biofuel can be extracted from fats or 
vegetable oils [149] and also algae [150]. The chemical kinetics 
development for biofuels has advanced to the next level by 
blending the n-heptane and methyl esters to represent rapeseed 
and soybean derived biodiesels, such as those developed by 
Herbinet et al. [148]. They reported that the ignition delay is in 
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Number of species/ Temperature Pressure Equivalence 
number of reactions range (K) range (bar) ratio range 
565/2540 550-1700 1-42 0.3-1.5 
159/770 625-1667 1-13.5 - 
57/290 600-1300 6-42 0.5-3.0 
26/52 700-1100 40-50 0.2-1.0 
71/360 - - - 
101/479 600-1400 10-40 0.22-2 
857/3606 550-1700 1-45 0.3-1.5 
45/69 - - 0.2-0.7 
32/55 750-900 30-50 0.2-0.6 
119/621 600-1300 20-50 1 
1389/5935 650-1200 3-50 0.5-1.3 
312/2469 650-1200 3-50 0.5-1.3 
48/67 700-1200 40 0.5-2.0 
53/325 1000-2500 0.01-10 0.1-5 
37/351 - - - 
55/278 900-1400 16-40 1 
19/15 1000-2000 1-30 0.5-1.5 
14/9 300-1800 40-60 0.7-1 
37/354 300-1800 0.02-100 0.5-3.0 
9/20 250-2500 0.05-87 0.2-6 
10/18 298-2700 0.05-87 0.2-6 
9/14 900-1400 1 0.5-5 
11/19 1000 1-30 0.5-4.5 
57/383 300-2500 1 0.5-2.0 
52/326 - - - 
112/484 - - - 
295/1498 800-1350 1-13 0.35-1.5 
53/156 650-1350 40-60 0.4-1.5 
3036/8555 800-1400 10 0.5 
115/460 700-1800 1-100 0.5-2 
426/2335 720-1700 0.04-80 0.6-1.7 
69/192 


good agreement with experiment using rapeseed oil. The mechan- 
ism can also be used for modeling the biofuel from various origins 
by adjusting the mole or mass fractions of the components 
(methyl decenoate, methyl-5-decenoate, methyl-9-decenoate 
and n-heptane) according to the actual biofuel blends. For more 
detail on chemical kinetic reaction mechanisms available for 
biofuels, one can refer to a review paper by Komninos and 
Rakopulous [47]. 


5.2. Numerical and CFD environments 


Numerical modeling can be categorized by the number of 
dimensions considered: zero-dimensional, quasi-dimensional 
and multi-dimensional. CFD is a multi-dimensional method 
whereby it resolves very small-scale zones and typically produces 
a more accurate result. A zero-dimensional model is the simplest 
model where there is only one independent variable, typically 
either time or, equivalently, crank angle, and an empirical heat 
release model is normally applied. The advantage of zero- 
dimensional model is that the run time is fast (computationally 
effective) and therefore it can be utilized for immediate use of 
engine design and analysis. However, the disadvantage of this 
method is that the simulated combustion duration is shorter than 
the actual duration due to inhomogeneities in reality [21]. The 
calculated results of zero-dimensional models produce rapid 
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pressure increases and very high heat release rates compared to 
actual HCCI engine combustion [143]. A quasi-dimensional model 
uses a turbulent sub-model for turbulent combustion and to 
derive a heat release model. Fig. 11 shows three different inter- 
acting regions for a quasi-dimensional model which is typically 
used to improve upon the zero-dimensional model. However, UHC 
emissions are over-predicted by 5-15% and CO emissions exhibit a 
50% error. This is due to the inability of the model to capture small 
temperature differences in the crevices [101,164]. 


5.3. Single-zone and multi-zone models 


A single-zone model is where the combustion chamber area is 
treated as one homogeneous block. A multi-zone model separates 
the combustion chamber into several zones and the zone distribu- 
tion depends on the type of the engine, either based upon flame 
propagation or homogeneous mixtures, thereby representing the 
inhomogeneity in the cylinder prior to combustion [21]. The zone 
distribution for HCCI engines is shown in Fig. 12(a) while Fig. 12 
(b) is a multi-zone model for SI engines. Multi-zone modeling of 
HCCI engines is organized in such a way because the heat 
generated in an HCCI engine comes from the core of the combus- 
tion chamber. This is different to SI engines where the zone is 
started from the spark plug (the location of the heat source), and 
propagated according to the flame front motion. 


Core Gas Unburned Boundary Layer 


Crevice Volume 


Fig. 11. General layout of the quasi-dimensional simulation showing the interacting 
adiabatic core, thermal boundary layer and crevice regions, reproduced from [101]. 


The single zone model has some limitations due to the 
assumption that the whole combustion chamber is treated as 
homogeneous. Peak cylinder pressure and rate of pressure rise can 
be over-predicted, as shown in Fig. 13(a). It also predicts a short 
burn duration and cannot accurately predict CO and HC emissions, 
which primarily depend on crevices [100]. Crevices and the 
boundary layer are the cold areas for HC and CO to react during 
combustion. On the other hand, the multi-zone model predicts the 
pressure trace and the peak cylinder pressure very well, as shown 
in Fig. 13(a), but it also cannot consider boundary layer effect and 
crevices, thus cannot predict CO and HC emissions [165]. Better 
agreement in combustion phasing is achieved when the single- 
zone model is coupled with detailed chemical kinetics and 
turbulent mixing effects [4], as shown in Fig. 11(b). 

If there is any inhomogeneity in the mixture, turbulence has an 
effect on combustion, while for a completely homogeneous 
mixture, turbulent mixing has little effect on the combustion heat 
release rate. However, it is suspected that there will be inhomo- 
geneities in the mixture as there is often insufficient time to mix 
down to the smallest relevant scales, which would cause turbulent 
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Fig. 13 (. a) Comparison between single-zone and multi-zone models with 
experiment [165]. (b) Comparison between experiments, single-zone model 
coupled with chemical kinetics only, and coupled with chemical kinetics and 
turbulent effects [4]. 


Fig. 12. Multi-zone model geometric configuration difference (a) HCCI engine, reproduced from [44] and (b) SI engines, reproduced from [166]. 
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mixing to directly affect reaction rates [4,6]. In actual combustion 
chambers, transportation of chemical species, heat transfer and 
heterogeneity of temperature and species concentration exist [43]. 
A multi-zone model of HCCI combustion predicted maximum 
pressure, burn duration, indicated efficiency and combustion 
efficiency with the worst error of 10% [165]. However, UHC and 
CO emissions were under-predicted; Aceves et al. [165] suggested 
that this might be due to some crevices not being considered in 
the analysis. They used a multi-zone model including zones in 
crevice regions producing more accurate results. 


5.4. Multi-dimensional models: CFD 


A multi-dimensional model solves the equations for mass, 
momentum, energy and species conservation to obtain more 
accurate results at the expense of computational time and 
resources. When solving detailed chemical kinetics, a sequential 
operation is used to reduce computational time by solving the 
flow field then the chemistry rather than attempting to satisfy 
both simultaneously. Two common approaches to multi- 
dimensional modeling are multi-zone models and CFD. A multi- 
zone model requires substantially less computational time and 
resources than CFD at the expense of accuracy. However, Aceves 
et al. [167] showed that a 40-zone model can successfully predict 
the effect of crevice geometry on HCCI combustion with combus- 
tion efficiency predicted within 5% error compared to experiment. 

A study has been performed using two-step processes in CFD to 
analyze combustion [10]. First, KIVA CFD was used for the effect of 
turbulence to solve the transport of all variables. Then the result 
from KIVA was used in the hydrodynamics, chemistry and trans- 
port (HCT) code to calculate the combustion parameters. This two- 
step method made it possible to obtain accurate predictions for 
turbulent combustion within a reasonable computational time 
[10]. The proposed approach using the CFD method can provide an 
accurate prediction of the combustion process and accounts for 
mixture inhomogeneities in both temperature and composition 
[168]. 


5.5. Effect of operating parameters on ignition delay time 


As discussed in Section 4, the ignition in HCCI engines can be 
controlled using a number of parameters. This can also be 
discussed using a numerical approach, where the chemical 
kinetics was used to investigate the ignition delay time for 
different fuels. Tanaka et al. [81] developed a reduced chemical 
kinetic model for a primary reference fuel (PRF), which is a blend 
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Fig. 14. The effect of initial temperature on ignition delay time for different fuels [81]. 
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Fig. 15. The effect of equivalence ratio on the ignition delay and burn rate for n- 
heptane, iso-octane and primary reference fuel (PRF) [81]. 
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Fig. 16. The effect of different pressures on the auto-ignition delay time for a 
stoichiometric mixture of n-heptane/air [170]. 


of n-heptane and iso-octane. They validated the model in a rapid 
compression machine (RCM) to represent HCCI combustion. The 
comparison was made between iso-octane, PRFs and n-heptane. 
The results in Figs. 14 and 15 show the ignition delay time for the 
selected fuel mechanisms with varying initial temperature and 
equivalence ratio. The results indicate that n-heptane has the 
fastest ignition delay time compared to PRFs and iso-octane. When 
the inlet temperature increases, the ignition delay time reduces 
further, where all fuels follow the same trend, as expected. The 
ignition delay time is also decreased when the equivalence ratio is 
increased towards the rich mixture, where the n-heptane has the 
fastest ignition delay time irrespective of the equivalence ratio 
compared to other fuels, as shown in Fig. 15. The n-heptane also 
has the highest burn rate, which shows its characteristic of 
susceptibility to auto-ignite. 

The results of ignition delay time between n-heptane and iso- 
octane by Tanaka et al. [81] is consistent with experiment. As 
discussed in Section 5.1, the n-heptane can easily be ignited using 
an inlet temperature of 313 K on a low CR engine. Gasoline, on the 
other hand, requires an inlet temperature of 453 K to auto-ignite 
for a high CR engine (CR=14.5) [169]. This shows that gasoline 
needs higher chemical activation energy compared with n-hep- 
tane. Different fuels require different inlet temperature to auto- 
ignite, as shown in Fig. 1, due to different ignition delay times. The 
inlet temperature requirement reduces as the compression ratio of 
the engine increases. This is also consistent with Olsson et al. [36], 
where they studied the effect of compression ratio on a natural gas 
fueled HCCI engine. They found that as the compression ratio 
increased, the required inlet temperature decreased, which 
reduced the ignition delay time. 

The effect of different inlet pressure using chemical kinetics is 
investigated by Puduppakkam et al. [170]. They studied the 
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ignition delay time between diesel and gasoline using detailed 
chemical reaction mechanisms. The n-heptane mechanism used as 
a surrogate fuel for diesel consists of 3809 species and 15,678 
reactions. The surrogate fuel for gasoline consists of 1833 species 
and 8764 reactions. Both surrogate fuels use a blend of several 
chemical classes, such as normal-, cyclo- and iso-alkanes, alkenes 
and aromatics. The result is shown in Fig. 16 for n-heptane, where 
the ignition delay time reduced as the inlet pressure increased. 
The pressurized intake air is able to increase the charge density 
and engine performance, and also can be used as an effective tool 
to control the ignition in HCCI engines, as discussed in Section 5.2. 

Fig. 17 shows the effect on ignition delay time when EGR is 
employed. It shows that the ignition delay time is slightly 
increased when EGR increases (reduction in O mole fraction). 
The result in Fig. 17 also shows that the ignition delay time reduces 
when the inlet pressure increases. The effect of EGR is also 
discussed by Christensen and Johansson [3], where they studied 
the mixture quality using iso-octane, ethanol and commercial- 
grade natural gas in a high CR HCCI engine. They reported that the 
ignition delay time also increased when EGR was used. This in 
return gives retarded combustion and slows the combustion rate, 
which is useful to reduce the combustion noise and also knocking. 

Knocking, as discussed in Section 2.8, is due to the instability in 
combustion phasing. Unstable combustion may cause knocking or 
misfiring (advanced or retarded ignition timing). EGR is able to 
retard the ignition timing and also to obtain higher loads, but too 
much retarded combustion causes misfiring especially for the fuel 
with longer burn duration [171]. The combustion phasing in HCCI 
engines is difficult to control at high loads because the stability 
limits between knocking and misfiring is narrow [172]. Maurya 
and Agarwal [173] reported that the unstable combustion is due to 
high cycle-to-cycle variation and to overcome this, they used 
closed loop control to monitor the combustion phasing. 

Olsson et al. [36] in their study of compression ratio influence 
on the maximum load of a HCCI engine reported that the 
combustion-phasing stability is due to the coupling between 
chemical kinetics and thermal problems. The gas temperature 
during combustion is influenced by the wall temperature or heat 
transfer rate, which affects the ignition timing. The instability in 
combustion becomes more severe at high loads, where the 
unstable combustion causes advanced or retarded combustion, 
as shown in Fig. 18. They employed hydrogen addition to control 
the ignition timing. The mass fraction of hydrogen was increased 
in order to advance the ignition. Once the wall temperature was 
substantially increased, the hydrogen amount was then reduced to 
keep the stability limit under control. Thus, a closed-loop control 
for the ignition timing was achieved by monitoring the wall 
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temperature (heat transfer rate) and the rate of hydrogen addition, 
in addition to the inlet temperature and fuel composition. 


5.6. Heat transfer model 


The heat transfer model plays an important role in combustion: 
it affects ignition timing; combustion duration; formation of CO, 
UHC and NO,; and the rate of in-cylinder pressure and tempera- 
ture rise, which is correlated to knock [174,175]. The influence of 
the heat loss through the combustion walls could be added in any 
of the numerical models discussed earlier. Some of the most 
commonly used heat transfer models are those devised by Anand, 
Woschni, Prandtl, Assanis and Hohenberg as in [176-178] and 
readers are advised to look at the work by Soyhan et al. [178], 
which explains the heat transfer equations by Assanis, Woschni 
and Hohenberg. From the study by Soyhan et al. [178], the 
Hohenberg model performs best in comparison with the Woschni 
and Assanis models: it predicted the in-cylinder pressure closest 
to the experiment, as shown in Fig. 19. 


5.7. Turbulence models: RANS vs. LES 
Turbulent flow occurs at high Reynolds numbers and, according to 


Tennekes and Lumley [179], the nature of the turbulence can be 
characterized as “irregularity, diffusivity, three-dimensional vorticity 
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Fig. 19. Heat transfer model comparison with experiment: (a) heat transfer coefficient and (b) in-cylinder temperature difference, reproduced from [178]. 
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fluctuations, large Reynolds number and continuum phenomenon 
(continuous in space and time)”. Currently there is no exact solution 
for most turbulent flows and this has led to the creation of turbulence 
models, where the useful methods to analyze the equation of motion 
is using statistical techniques [180,181]. The simplest statistical tech- 
nique is to consider the average of flow variables over time. In HCCI 
engines, a piston-crown design is important to create turbulent 
behavior inside the combustion chamber. A study by Kong et al. [64] 
found that the square-bowl piston generates higher turbulence levels 
and results in a longer combustion duration, where the turbulence 
affects the combustion through its influence on the wall heat transfer 
and property transport. This is supported by Aceves et al. [182] with 
their experimental data, where the square-bowl piston resulted in a 
longer combustion duration but with higher UHC and CO emissions 
and lower combustion efficiency as a result of the thicker boundary 
layer. Therefore, one has to pay particular attention when designing 
the piston crown to achieve the desired combustion behavior. Another 
study by Kong and Reitz [4] revealed that including the turbulent 
mixing effects gave better predictions for in-cylinder pressure com- 
pared to the chemistry alone, as shown in Fig. 20. The link between 


chemistry and mixing effects with turbulence models is illustrated in 
Fig. 21. 

Reynolds-Averaged Navier-Stokes (RANS) turbulence models 
are used in transient fluid flows, whereby the flow variables are 
decomposed into average and fluctuating quantities over regions 
in physical space that are of the order of the integral length scale 
[180]. Different classes of RANS models exist: zero-equation 
models, one-equation models, two-equation models and stress- 
equation models [185]. RANS models have their disadvantages, 
whereby: (1) they are not able to properly predict simple free- 
shear layers and are known to be inadequate in complex flows 
[185,186]; (2) they are unable to capture dynamically evolving 
fine-scale vortical and scalar structures due to turbulence model 
induced dissipation on the in-cylinder unsteady mixing process 
[187]; (3) the k-e model has poor performance when used in 
unsteady flows, where the model over-predicts the turbulent 
length scales resulting in a high wall shear stress [188,189]. This 
is supported by Celik et al. [190], where it performs poorly when 
applied to IC engines with the predictions being less accurate than 
desired. 

Large Eddy Simulation (LES), on the other hand, is a turbulence 
model that computes more turbulent length scales than RANS, 
with the turbulent energy within the inertial sub-range separated 
into resolved large-scale eddies and unresolved small-scale eddies. 
Resolved large-scale eddies are computed directly by the discre- 
tized equations, while small-scale unresolved eddies, which con- 
tain only a small fraction of the turbulent kinetic energy, are 
modeled through subgrid scale models. LES equations are obtained 
by appropriate filtering of the Navier-Stokes equations and the 
energy equations over a three-dimensional space [190]. In com- 
plex flows for nonreactive and reactive system studies, LES 
predicts the scalar mixing process and dissipation rates with 
considerably improved accuracy compared to RANS [191]. It offers 
a more realistic representation of the in-cylinder turbulent flow 
and delivers accurate predictions for non-premixed combustion 
[184,192]. 

However, LES has difficulties in the near-wall region, where the 
mesh should be fine enough to compute a correct shear stress 
[184,193]. Furthermore, results from LES may introduce serious 
errors into simulations due to both aliasing and truncation errors 
of low-order schemes, which can degrade LES computations | 180]. 
Therefore, LES requires greater mesh refinement compared to 
RANS, which increases the computational cost. However, mesh 
or time-step refinement in RANS improves the numerical accuracy, 
but does not improve the scales that are resolved. Instead, it 
converges to an exact solution of the filtered Navier-Stokes 
equations, while LES converges towards the unfiltered Navier- 
Stokes equations [192]. 
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Fig. 21. Illustration of an overall model for transient simulations [180,183,184]. 


A hybrid approach between RANS and LES has been created in 
order to resolve the disadvantages of these two models. Alfonsi [185] 
has categorized this approach into zonal decomposition, nonlinear 
disturbance equations and universal modeling. Benarafa et al. [193] in 
their study on the hybrid approach, which they called RANS/LES 
coupling, shows that this approach seems to correctly force the 
averaged LES velocity to reach the correct steady RANS velocity field 
even at high Reynolds numbers. The fluctuations are globally 
improved despite the use of very coarse meshes. 

Another approach which is similar to the hybrid RANS-LES 
model is called Detached Eddy Simulation (DES). DES, which is 
already included in commercial CFD packages, implements the 
RANS approach near the wall region, while using the LES techni- 
que in the outer detached eddies. The main objective of using this 
technique is to avoid using high resolution in the near-wall region, 
and therefore, reduce computational cost. Instead, algebraic 
boundary-layer models (which have proven to be very successful 
in the RANS context) are implemented. 


5.8. Mixing models for non-premixed combustion 
Turbulent mixing will take place when there is inhomogeneity 


in the mixture and depends on the geometry of the piston crown. 
Fuels and oxidizers are required to be mixed at molecular levels to 


initiate combustion, while the molecular mixing will take place at 
the interface between small eddies [180,191]. This process will 
heavily involve the chemical reactions, which leads to the diffusion 
process. 

The Euclidean Minimum Spanning Tree (EMST) [194] mixing 
model is designed to overcome shortcomings of simpler turbulent 
mixing models, and was successful in diffusion flame tests 
[194,195]. The problem with simpler models is that the mixing 
occurs between particles across the reaction zone, which is not 
local. The EMST model only allows the particles to mix between 
their immediate neighbors in composition space. The model was 
constructed by connecting particles to form branches of trees over 
the timestep. The particles enter or leave the cell according to their 
velocity. Each of the particles has its own age property and this 
determines whether the particle belongs to the mixing or non- 
mixing state. If the particle age property is positive, it belongs to 
the mixing state, while it is not if the age property is zero or 
negative. The EMST is formed on the mixing state and changes 
discontinuously in time due to the particles entering or leaving the 
mixing state. 

Mitarai et al. [195] conducted a performance comparison 
between EMST, Interaction by Exchange with the Mean (IEM) 
[196] and the Modified Curl’s (MC) [197] mixing models in RANS 
and LES environments and found the EMST mixing model yielded 
significantly better results than the IEM and MC mixing models. 
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Fig. 22. Scatter plots of temperature vs. mixture fraction from experiment and simulations using LES: (a) experimental data, (b) IEM, (c) MC and (d) EMST. All figures 


reproduced from [198]. 


It was found that the EMST mixing model performs better in an 
LES environment than in a RANS environment. This is because of 
the ability of LES to resolve finer scales compared to RANS and 
therefore particle interactions are over finer scales [198]. The 
EMST mixing model has also successfully predicted the appear- 
ance and disappearance of cold particles in LES, where it has failed 
in RANS. A study of turbulence and chemistry interaction per- 
formed by Bisetti and Chen [198] showed that the EMST in the LES 
environment performed better than the IEM and MC in predicting 
the temperature conditioned on the mixture fraction, as shown in 
Fig. 22. 

Bisetti et al. [199] studied turbulence and chemistry interac- 
tions in HCCI engines using the same models as Mitarai et al. [195]. 
They found that the EMST mixing model performed better than 
the other mixing models. They noted that the LES simulation gives 
a promising step in HCCI engine modeling under high levels of 
stratification. 

The Probability Density Function (PDF) [200] method offers a 
good advantage in modeling turbulent reactive flows, where it has 
the ability to capture strong turbulence-chemistry interactions in 
the mixing models, such as IEM, MC or EMST [201]. But, this model 
is hard to solve because it has many dimensions: one for each 
chemical species in a particular set of chemical reactions. There- 
fore, it will end up with n dimensions for n chemical species and if 
there are m nodes per dimension, the total number of nodes will 
be m”, resulting in an enormous computational cost, which is why 
it is always solved stochastically. The Conditional Moment Closure 
(CMC) [202] model on the other hand, was designed to handle 
turbulence chemistry interactions, where large chemical mechan- 
isms can be used at modest computational cost. The conditional 


fluctuations of the reactive components are to be neglected when 
calculating the chemical source term for first-order CMC: the value 
is significantly smaller than the unconditional fluctuations [203]. 
Bushe and Steiner [204] implemented the CMC model in the LES 
environment for non-premixed turbulent reacting flows. They 
noted that the resolution constraints were reduced because the 
chemical reactions are resolved in the mixing space. However, the 
CMC model has a disadvantage in modeling the properties’ 
fluctuations since it has only one conditioning variable [205]. 

The MMC model [205] was designed to solve the difficulties faced 
by the PDF method and the CMC model. It combines the precepts of 
PDF modeling (in modeling the major species) and CMC (for the 
minor species) while treating the complete set of dependent variables 
equally [205]. The major species are solved by mapping them to a set 
of prescribed reference variables (often standard Gaussian in the 
conventional description), while fluctuations of minor species are 
either ignored (conditional MMC) or treated with conventional mixing 
models (probabilistic MMC) [206]. 

An investigation of MMC, studying the same case as Mitarai 
et al. [195], reported that the model is capable of predicting the 
mean temperature rise with reasonable fidelity and follows the 
temperature history at stoichiometric mixture fraction better than 
the other models [207], as shown in Fig. 23. A study by Cleary and 
Kronenburg [206] for the MMC model in turbulent diffusion 
flames shows an improved accuracy in local extinction of chemical 
time scales compared to the CMC model. Vogiatzaki et al. [208] 
tested the MMC model for inhomogeneous reactive flows in 
turbulent jet diffusion flames. The model was constructed using 
a one-dimensional reference space with mixture fraction used as a 
basis for the reference variable. They found that the temperature 
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Fig. 23. Temperature at stoichiometric mixture fraction for various mixing models 
[195,207]. 


was predicted well compared to the experiment at all locations. In 
addition, Vogiatzaki et al. [209] mentioned that the MMC model 
could be easily included in an LES turbulent model. This was done 
by Cleary et al. [210] where they used the MMC model with an LES 
scheme for turbulent diffusion flames. The LES model was used to 
solve for the turbulent velocity field and the reference mixture 
fraction, while the stochastic MMC model solved the reactive 
scalars. They concluded that the conditional averages of tempera- 
ture and species were in good agreement with the experiment. 

The advantages of using the MMC model [209,211] are: (1) the 
simplicity of modeling the MMC diffusion coefficients against other 
models, (2) conditional scalar dissipation rates appears in closed form 
and (3) no need to presume the mixture fraction PDF because it is 
computed from the modeling of the relationship between mixture 
fraction and reference variable. However, there is difficulty in model- 
ing the MMC, where the model cannot generate fluctuations around 
the conditional mean for the joint PDF, and an additional transport 
equation is necessary to impose the conditional fluctuations [212]. In 
other words, the MMC method looks to be a promising step in 
developing a mixing model for HCCI engines. 


6. Conclusion 


It can be clearly seen that the hydrogen-diesel combination in 
HCCI engines shows better results for most major quantities, while 
producing a higher efficiency compared to the single diesel mode and 
the natural gas—diesel mode. A higher BTE will help reduce the fuel 
consumption of the engine and it might be able to match the fuel 
consumption of a hybrid engine. If the HCCI engine is used in a hybrid 
configuration, it might help reduce the fuel consumption even further. 
It looks like a promising choice of engines in the near future. 

The HCCI engine has low emissions levels of NO,, soot and 
particulates. However, HCCI engines still have unresolved issues, 
which are knocking and high levels of unburned HC and CO 
emissions. If hydrogen is added to diesel-air mixture, knocking 
will take place once the energy ratio is more than 16%. It is 
expected that the knocking in HCCI engines is influenced by 
chemical kinetics, where the combustion is no longer controlled 
by the injection or the spark timing. Therefore, it is hypothesized 
that the knocking in HCCI engines could be due to the region with 
the highest OH concentrations in the high temperature zone, 
where the highest temperature region will initiate instantaneous 
local combustion resulting in a local high heat release rate. Further 
studies have to be performed in order to solve these remaining 
issues. To achieve this, the numerical method shows a great 
advantage over experiments to simulate the combustion behavior 
in terms of cost and time. To this end, a simulation model has to be 


developed to investigate the behavior and, once complete, it has to 
be validated against experiments. 

The multi-zone numerical method combined with advanced 
turbulent mixing models shows promising results, compared with 
the single-zone model. CFD offers greater accuracy than the multi- 
zone model at the expense of computational cost. But, one could 
easily implement the LES or DES model with CFD. The DES model 
will be useful in the near future, where it combines the LES and 
RANS models together with improved accuracy using the former 
where needed and improved computational efforts using the latter 
where possible. A mixing model has been developed to study the 
mixing behavior at molecular levels. The performance of EMST, MC 
and IEM models were compared and it showed that the EMST 
model performed best. A study of the CMC model showed that it 
has the ability to capture all phases of combustion, i.e. ignition, 
flame propagation and a diffusion flame [213]. A CMC model in an 
LES environment has been investigated, and it shows that the 
model has the capability to predict extinction and ignition phe- 
nomena [204]. The MMC model is created by combining the 
advantages of the CMC and PDF models. Therefore, the MMC 
model is expected to deliver a better result if it is implemented in a 
DES environment, as both models exist to use the advantage of 
other models. 

There are several methods to achieve HCCI configurations through 
experiments. Using a direct injection system seems to require few 
modifications, but care is required with the associated disadvantages. 
If port injection is to be used, one might have to install the atomizer 
and heater in the inlet port, which leads to additional costs. The main 
objective of the HCCI engine is to achieve a high compression ratio, 
which provides similar efficiency to the CI engine or better, with 
emissions levels at least as good as SI engines. Therefore, an optimized 
experimental method would be able to solve the HCCI engine 
problems, with the help of simulations. The major problem yet to 
be solved is knocking, thus it is imperative to develop a model so that 
HCCI engines can be used in the near future. 
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